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Abstract 

We compute one-loop corrections to the annihilation of non-relativistic particles 
X due to the exchange of a (gauge or Higgs) boson 99 with mass /i in the initial state. 
In the limit ^ this leads to the "Sommerfeld enhancement" of the annihilation 
cross section. However, here we are interested in the case /x < m^,, where the one-loop 
corrections are well-behaved, but can still be sizable. We find simple and accurate 
expressions for annihilation from both S— and P— wave initial states; they differ from 
each other if 7^ 0. In order to apply our results to the calculation of the relic density 
of Weakly Interacting Massive Particles (WIMPs), we describe how to compute the 
thermal average of the corrected cross sections. We apply this formalism to scalar 
and Dirac fermion singlet WIMPs, and show that the corrections are always very 
small in the former case, but can be very large in the latter. Moreover, in the context 
of the Minimal Supersymmetric Standard Model, these corrections can decrease the 
relic density of neutralinos by more than 1%, if the lightest neutralino is a strongly 
mixed state. 
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1 Introduction 



The existence of non-baryonic Dark Matter (DM) in the universe is by now well established 
[1]. Recent observations determine the universal average DM density quite accurately. 
The exact value and its uncertainty depend somewhat on the assumptions made in the 
fit ("priors"). For example, within a minimal ACDM model, and combining data on the 
cosmic microwave background (CMB) anisotropies with observations of supernovae of type 
la and with analyses of baryon acoustic oscillations, one finds [2]: 

(]cDM/i^ = 0.1131 ±0.0034. (1) 

Here I^cdm is the energy density of cold dark matter in units of the critical density, and 
h is the scaled Hubble parameter such that Hq = lOO/ikmsec"^ Mpc~^ where Hq is the 
current Hubble parameter. Introducing additional parameters in the fit can increase the 
allowed range; for example, significantly larger values of QcBuh'^ are allowed if the number 
of particles that were relativistic when the CMB decoupled is kept free |2]. However, the 
minimal model describes the data well. Moreover, quite soon data from the Planck satellite 
are expected to reduce the error on flcDuh'^ to the level of 1.5% using CMB measurements 
alone [3]. 

From the particle physics point of view. Weakly Interacting Massive Particle (WIMPs) 
are among the most attractive DM candidates. In standard cosmology their thermal relic 
density is naturally of the right order of magnitude. Owing to their weak interactions, they 
can be probed through both direct and indirect detection experiments [Tj, although no 
convincing signal has yet been found. Moreover, the existence of WIMPs is independently 
motivated in many extensions of the Standard Model (SM) of particle physics, the most 
prominent example being extensions based on supersymmetry (SUSY). 

In order to fully exploit the precision of present cosmological data, the theoretical error 
on the prediction for the WIMP relic density should not exceed the error inferred from ob- 
servations. For a given cosmological model, the WIMP relic density is determined uniquely 
by its annihilation cross section into SM particles, if all WIMPs were produced thermally. 
This requires that the post-inflationary Universe was sufficiently hot, with temperature 
exceeding about 5% of the WIMP mass. Moreover, one has to assume that no other, heav- 
ier particles decayed into WIMPs after WIMP decoupling. Under these conditions, which 
are satisfied for standard cosmology, the uncertainty of the current WIMP relic density is 
essentially given by the uncertainty of the WIMP annihilation cross section. In order to 
calculate this cross section to percent level accuracy, at least leading radiative corrections 
will have to be included]^ 

In this paper we calculate one class of potentially sizable radiative corrections, which 
are due to the exchange of a boson between the WIMPs prior to their annihilation. If the 
mass fi of the exchanged boson vanishes, the one-loop expression diverges in the limit of 
vanishing relative velocity v between the annihilating WIMPs. These large "Sommerfeld" 

'''These corrections also contribute to the cross sections for WIMP annihilation into SM particles in our 
galaxy at present times, which affect the size of indirect WIMP detection signals. However, currently there 
are very large uncertainties in the backgrounds to these signals; in case of charged particles, propagation 
through the galaxy adds additional uncertainty. Percent level correction to signals for indirect WIMP 
detection are therefore not significant. 
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corrections then have to be re-summed to all orders in the relevant coupling [U [3, [6] (for 
earlier, related work see [7]). However, since WIMPs couple neither to photons nor to 
gluons, a new light boson with sizable coupling to the WIMPs, but not to SM particles, 
has to be introduced. The required hierarchy fi <^ m^, where is the mass of the WIMP, 
then raises naturalness issues. 

Here we instead study the case fi < m^. Examples are the exchange of the hght Higgs 
boson in the minimal supersymmetric extension of the SM (MSSM), and the exchange of 
the Z boson in most "generic" WIMP models. If /i < and f ^ 1 these corrections can 
be significantly bigger than the "generic" estimate a/vr ~ (a few) xlO~^, while remaining 
small enough to permit simple one-loop calculations (and without threatening unitarity 
[Sj). We calculate these corrections in the non-relativistic limit, following lengo [5]. Since 
WIMPs decouple at temperature T ~ an expansion in the relative velocity usually 

(but not always [5]) converges quite fast, and can therefore be used in the calculation of 
the one-loop correction. We also show how to compute the relevant thermal average over 
the corrected cross sections. We first apply our formalism to simple models with scalar or 
fermionic SM singlet WIMPs; the corrections are always small in the former case, but can 
be sizable in the latter scenario. We then analyze neutralino annihilation in the MSSM, 
and show that these corrections can reduce the relic density of strongly mixed neutralinos 
by more than 1%, comparable to the projected uncertainty of the value to be inferred from 
Planck data. 

The rest of this paper is organized as follows. In Sec. [2] we introduce the formalism, 
which was initially proposed to treat non-perturbative Sommerfeld enhancement [Ej. In 
particular, we factorize the correction to the amplitude due to the boson exchange between 
two incoming particles. In Sec. [3] we briefly review the standard calculation of dark matter 
relic density, and describe the thermal averaging of the corrected cross section. Numerical 
results for three WIMP models will be presented in Sec. |H Finally we will summarize. 
Some details of our numerical procedure are described in Appendix A, while Appendix B 
shows that our approximate treatment indeed reproduces the leading terms of an exact 
calculation of radiative corrections associated with the initial state in a simple scalar model. 

2 Correction to the annihilation amphtudes 

Consider the annihilation of two WIMPs x iiito two SM particles: 

xiPi) + xiP2)^X^{p[)+XM). (2) 

Generic tree-level diagrams contributing to this process have the form shown in Fig. [H 
We want to compute one-loop corrections of the kind shown in Fig. [21 where a boson (f is 
exchanged between the WIMPs before they annihilate, by adapting the formalism of lengo 
[5], [To]. We assume that x is a Majorana fermion; however, in the non-relativistic limit this 
will be relevant only for the case where the exchanged boson has axial vector couplings 
(see below). 

In this formalism (f exchange and x annihilation are factorized; the former can then 
also be understood as re-scattering of the incoming WIMPs prior to their annihilation. 
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This factorization [ip exchange before x annihilation) is only expected to work in the non- 
relativistic limit. Moreover, it requires the virtuality of the ip propagator to be (much) 
smaller than that of the particle exchanged in x annihilation. The latter can always be 
satisfied for /i ^ m^, where // is the mass of the exchanged boson, but it can also be 
satisfied for yU ~ if the WIMPs annihilate through the exchange of a particle Y with 
mass My ^ m^^. However, we will see that the corrections become small if /i > m^. 
Note also that these corrections do not capture UV effects like the renormalization of the 
tree-level couplings of x- 




Figure 1: Tree-level WIMP annihila- 
tion. The grey blob represents the ex- 
change of some particle in the s—,t— 
or M— channel. 



Figure 2: The one-loop corrections 
we are considering in this article. The 
meaning of the grey blob is the same 
as in Fig. 1. 



Let P = {pi+ p2)/2 = {p[ + ^'2)72, p = {pi- p2)/2 and p' = {p[ - p'.^jl; recall that 
Pi, P2 and p'y^ p^ are the four-momenta in the initial and final state, respectively. In the 
cms we have Pq - 



fp' -\- m^, P = and po = 0. The initial state kinematics is thus fixed 
by p. We write the one-loop corrected amplitude for annihilation from the partial wave 
denoted by L as 

Al{\p\,p') = Ao,l{\p\,p') + 5Al{\p\,p') , (3) 

where Aq^l and 5Al denote the tree-level amplitude and the one-loop correction term, 
respectively. We are only interested in the cases L = S and L = P. The contribution of 
higher partial waves is strongly suppressed during WIMP freeze-out. 

The corrections can be calculated starting from the observation that annihilation from 
an S* — (P— )wave initial state can be described by a pseudoscalar (scalar) x current [TT] . 
This allows to write the correction as a standard vertex (three-point function) correction. 
Let us begin with the simple case of scalar boson exchange: 



sAlUpI^p') 



ig'^v{p2) 



(75)"^ 



(27r)4 {q - P)2 -ml + ie' {q + Pf - + ie 
1 



X 



{p-qf 



le 



-Ao^Li\q\,v')u{j>i) . 



(4) 



Here g is the strength of the coupling between the boson and the WIMP, = 1 (0) for 
annihilation from an S* — (P— )wave initial state, and the reduced tree-level amplitude 
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Aq^l describes the blob in Fig. 1 [except for the factor 75 in case of 5— wave annihilation, 
which appears explicitly in Eq.(jl])] as well as the final state|f| We will later show how to 
treat the exchange of spin-1 bosons. 

As a first simplification, one then uses the fact that the full relativistic boson propagator 
satisfies — /i^) = — l/(/c^ + /x^) [1 — A;q/(/c^ — yU.^)]; the second term can be omitted 

in the non-relativistic limit, where the energy exchange is much smaller than the three- 
momentum exchange. Moreover, to leading order in the non-relativistic expansion, the 
g— dependence of the reduced amplitude can be neglected, i.e. the factor Aq^l can be 
pulled out of the integralj§ We then perform the integrals in Eq. (jl]) starting from the go 
integration. We are looking for poles in the lower half-plane, with residues that diverge in 
the limit g, p — )■ [5] . This gives 

r ^ /N 2-f \ f {.i - P + ^y) {l^T'^ {i + P + '^y) ~A \ ,r\ 

) ^ 9 v{v2) / TTrk^- ' . 1 J 3^^o,l(p') u{vi) , (5) 



(2vr)3 8a;Po(w-Po) + 

where u = ^(p + m^; in the numerator one should take go = co — Pq. For small 3-momenta, 
we have 

uj-Po^ — — ; (6) 

note that this vanishes for p, g — )■ 0, as advertised. 

To zeroth order in the non-relativistic expansion we can set g — )■ in the numerator 
of Eq.([5]). We see that this gives a non-vanishing result only if = 1, i.e. if a 75 matrix is 
present; recall that this corresponds to annihilation from an S*— wave. To this order we can 
replace p = Pq'J^ acting to the right (on the m— spinor) by m^. The numerator of Eq.([5]) 
then reduces to A'-f^m'^^. Note that the factor of 75 is required in order to be able to access 
the large component of v{p2)] this is most easily seen in the Dirac representation. Moreover, 
the factor uPq in the denominator of Eq.([5]) can be replaced by m^, up to corrections which 
are of second order in three-momenta. 

To summarize, we have made three approximations: 

1. We ignored the energy dependence of the (p propagator. 

2. In the go integral we only kept the pole with the leading residue (in the non- 
relativistic limit). 

3. We ignore all g dependence in the numerator (for annihilation from the S"— wave). 



^Notc that Eq.Q is one-loop exact, if the q dependence of Aq l is kept. For a full non-perturbative 
treatment the complete reduced amplitude should appear again in the integral on the right-hand side 
[5]. Recall, however, that we are only interested in one-loop corrections, in which case we may use Aq^l 
in the integrand. 

®In the P— wave case, the nontrivial dependence on the initial state three-momentum stems from the 
spinors describing the initial state and the Dirac structures shown explicitly in Eq. Q , not from the reduced 
amplitude. Aq^l in fact does not depend on |(7| if w annihilation proceeds through an s— channel diagram. 
For t— or u— channel annihilation wc have to assume that the particle exchanged in the annihilation process 
is significantly more off-shell than |g| ~ \p\. 
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Note that these approximations have to be taken simuhaneously in order to get a UV- 
finite resuh. One may worry that this gives a rather poor approximation of the exact vertex 
correction even in cases where the latter are finite, as in a purely scalar theory. In Appendix 
B we show that this approximation can indeed differ substantially from the exact vertex 
correction if /i <^ m^. Note that the exact vertex correction by itself becomes IR-divergent 
for yU = 0. This leads to terms cx: log(m^//i) which appear in the exact vertex correction, 
but not in our approximation. However, these terms are canceled by real emission diagrams 
and wave function corrections, which have to be included in order to obtain an IR-finite 
result for /i = 0. In Appendix B we show that, at least for a simple scalar model, our 
approximation does accurately reproduce the exact radiative correction associated with 
the initial state, whenever these corrections are large. 

Let us therefore proceed with our calculation, which does not require additional ap- 
proximations. The angular integrations are straightforward. One is then left with a single 
integral to describe the correction to S*— wave annihilation: 

(27r)'^ \p\ Jo q^-p^ {\p\-\q\) +1^ 

Note that we have absorbed the spinors into the full tree-level amplitude Aq^s- The one- 
loop correction to the cross section emerges from the interference between the correction 
6As and the tree-level term ^0,5. We can thus write 

5^5(b1,pOli-ioop = £^lM\p\,pVs{r) , (8) 

where v is the relative velocity between the two annihilating WIMPs in their center of 
mass frame (i.e., |p| = m^v/2), and we have defined the function 



Isir) = 3fJe 



' X , (1 + x)^ + r 

T n 

x'^ — \ (1 — xV + r 



(9) 



Here x = \q\/\p\ and r = fi^/jf. 

So far we have assumed that a spin-0 boson with scalar coupling is exchanged between 
the annihilating WIMPs. In order to treat more general cases, we rewrite the numerator 
ofEq.(l5l)as 

M = Tii-p + m^){^,m + P + m^)T, (10) 

where F describes the Dirac structure of the (pxx coupling and F = 7°F"''7° its Dirac 
conjugate. Scalar exchange corresponds to F = F = 1. 

It is easy to see that pseudoscalar exchange, F = — F = 75, does not lead to enhanced 
contributions. For example, forn^, = 1 and g — one finds a result oc 75 {rriy. — /^)^, which 
is 0{f). 

Vector exchange corresponds to F = 7^, F = 7j,, where the Lorentz index u has to be 
summed. For = 1 and g — j- this gives 

A/'vcctor = -2m^-f5Y {rn^ + P)lu = -'^m^lb (Sm^ - f^) . 

Again replacing p by this leads to the same result as for scalar exchange, except for 
an overall sign. However, this sign is compensated by the extra minus sign in the vector 
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boson propagator. We thus reproduce the well-known result that in the non-relativistic 
limit, the exchange of a vector boson has the same effect as that of a scalar boson. 

Finally, axial vector exchange is described by setting T = 'j'^'j5, f = 7^75, where 
summation over u is again implied. In the S*— wave case, where = 1 and g — )■ in the 
numerator, this gives 

A4xial vector = 757^' iP " f^xf 1^ = 4"^x75 (2m^ + p) . 

Again replacing p by m^, and accounting for the minus sign in the spin-1 propagator, we 
find that axial vector exchange differs from scalar or vector exchange by a factor of —3. 

This does not seem to have been noticed in the recent literature. We therefore checked 
it in the limit where ip exchange can be treated as re-scattering of the incoming WIMPs, 
leading again to two on-shell WIMPs. To leading order in velocity expansion we are in- 
terested in the limit of vanishing momentum exchange. However, we have to keep in mind 
that the two WIMPs will have to annihilate through a 75 vertex (for the S"— wave case). 
This requires that the non-relativistic u and v spinors have the same spin. The rescat- 
tering process is then described by the quantity (we omit the bosonic propagator and all 
couplings) 

^res = ^V{p2, s)Tv{p2, S')u{pi, s')Tu{pi, S) , (11) 
s' 

where s,s' describe the spin. Scalar boson exchange again corresponds to P = P = 1. In 
this case only s' = s contributes in the non-relativistic limit, where (in the cms) pi = p2 — 
(m^, 0), and one has Ares, scalar = ~4m^. One immediately sees that pseudoscalar exchange 
only contributes at OijP'). 

In case of vector exchange, only P = 7° contributes to 0{p^) [5]. This then again 
requires s = s' in Eq. ffTT]) . giving A^^s, vector = 4m^. Remembering the additional minus sign 
in the spin-1 propagator we therefore again find that vector boson exchange contributes 
the same way as scalar exchange. 

Finally, for axial vector exchange, only P = 7^75 {i = 1,2,3) contribute to 0{p^). In 
this case Aj-es receives non-vanishing contributions both from s = s' and from s = —s', 
so that Aj-cs axial vector = — 12m^. Again including the minus sign from the propagator, we 

' A. 

reproduce our earlier result that axial vector exchange differs from scalar or vector exchange 
by a factor of —3. Note also that axial vector exchange describes an interaction between 
the spins of the two WIMPs. Such interactions are not suppressed at low velocities. 

Now let us discuss P— wave annihilation, which corresponds to = in Eq.([5]). In 
this case setting q = would lead to a result which is of 0{jf), i.e. of second order in 
the non-relativistic expansion. The leading term is the one linear in q; the numerator of 
Eq.([S]) then becomes —im^q-^. Note that the 7^ again allow to access the large component 
of v{p2). The proper form of the correction term can then most easily be obtained using 
trace techniques, by dividing the 1-loop correction term SApAqp by the tree-level result 
\Aop\ . We find that the P —wave correction term differs from the S*— wave term by a factor 
p- q/p^ inside the momentum integral. Performing the angular integrals, we can write this 
in the form of Eq.([H]), with a new function describing the correction for annihilation from 
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a P— wave initial state: 



Ip = 3fte 



oo 



dx 



X- 



:2 - 1 



2x2 



1 + 



X + 1 + r , (x + 1) + r 

m 

Ax [x — 1)^ + r 



1 



(12) 



where r = fi'^/p^ as above. 

In order to treat the exchange of other bosons, non-trivial Dirac structures F, T again 
have to be introduced in Eq.([5]), now for the case = 0. Proceeding as above, we find that 
pseudoscalar exchange does not lead to a large correction, while vector and axial vector 
exchange give the same correction as the exchange of a scalar boson. At first glance it 
may seem surprising that now axial vector exchange gives the same, positive, contribution; 
recall that in case of annihilation from an S*— wave, axial vector exchange differed by a 
factor of —3. The difference can be understood from the observation that axial vector 
exchange leads to a spin-spin interaction of the form 4si ■ S2, where si^2 are the spins of the 
two WIMPs. This can be evaluated using 2si ■ S2 = S"^ — si — S2, where 5* = si + 5*2 is the 
total spin. In case of Majorana WIMPs, annihilation from an S"— wave requires [H] S = 0, 
leading to 4si ■ S2 = —3. For P— wave annihilation, we need [H] 5" = 1, giving 4si ■ S2 = 1. 
Note the relative factor of —3 between these two results. 

At this point a comment on other WIMPs (than Majorana fermions) is in order. For a 
Dirac fermion-antifermion pair, there is no strict correspondence between the total spin and 
the orbital angular momentum. There will also be S*— wave states with S = 1 (as the J/ip 
family of quarkonia), as well as P— wave states with S* = 0. In this case the proper factor 
in front of the axial vector correction will depend on the spin state, i.e. it will no longer 
be completely process independent in this case. However, results for scalar, pseudoscalar 
and vector exchange are the same as for Majorana WIMPsjll Finally, there is no such thing 
as an axial vector coupling to scalars, but our results for vector exchange apply to scalar 
WIMPs as welS Scalar exchange now involves a trilinear scalar interaction /, which has 
dimension of mass. Our results can describe this situation as well, with g = //(2m^). 

The integrals in Eqs. Qj and f|T2l) should be understood as principal value integrals, in 
order to treat the pole at x = 1. In case r ^ 1 the integrals can be computed analytically, 
by expanding the logarithms in inverse powers of r. Moreover, in the limit r — )■ both Is 
and Ip approach 7r2/2, thereby reproducing the well-known result [12] that the one-loop 
"Sommerfeld factor" for massless boson exchange is the same for S— and P— partial waves. 
We also found accurate numerical expressions for small and moderate r. Altogether, the 
correction factors can be described by 



^Strictly speaking, Majorana fermions do not have diagonal vector couplings. They can, however, have 
vector couplings to other Majorana states. Our result is applicable to this situation in the limit where the 
mass of the second Majorana fermion approaches that of the annihilating WIMP, 
vector coupling exists only for a complex scalar, which can carry a charge. 




(13) 
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The approximations for large and small r intersect at r ~ 5.6 (4.2) for Is [Ip)] the two 
approximations for Ip intersect a second time at r ~ 6.3. 




Figure 3: Comparison between the exact Is (solid, black) and its approximations for large 
r (dotted, green) and small r (dashed, red). 

In Fig. 131 the exact function Is and its approximations are shown. By switching from 
the low— r to the high— r expression at the intersection point one reproduces the exact nu- 
merical result to better than 4% for all values of r. In case of Ip (not shown) the large— r 
approximation overshoots rather than undershoots the exact result for 10"'^ < r < 1. This 
leads to slightly larger discrepancies, of up to 6%, between the exact Ip and its approxi- 
mations at r ~ 5 where the two approximations intersect. For the purpose of calculating 
the WIMP relic density a relative error on the correction of 6% is quite acceptable. 

In Fig. |4]we show I sir) and Ip{r) as black and red (grey) lines, respectively. As noted 
earlier, the two functions coincide for massless exchange bosons, r = 0. For nonvanishing 
boson mass, r > 0, the S*— wave contribution is larger than the P— wave one, by up to 
a factor of 3 at large r; see Eqs. fll3p . Note that for any finite mass of the exchanged 
boson, fJ, ^ 0, the zero-velocity limit |p| — >■ corresponds to r — >■ cx). Eqs. lIT^ show that 
asymptotically Is,p{.'f cxd) oc l/i/r = Eq.® then shows that the corrections 

approach constant values of order g^m^j {4:Hfi) for |p| — )■ 0, if /i 7^ 0. Such corrections will 
threaten the convergence of perturbation theory only if the WIMP and boson masses differ 
by a loop factor; in the technically more natural case where the boson mass lies a factor 
of a few below that of the WIMP, we still find a significant enhancement relative to the 
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Figure 4: Comparison between Is describing the correction for annihilation from an 
S*— wave (black) and Ip describing the correction for annihilation from a P— wave (red 
or grey). 

naive expectation that corrections should be of order g"^ / {IGn'^) . We finally note that our 
numerical results are consistent with those in Ref. [6]. 

3 Dark Matter relic density 

In this section we describe how to calculate the loop-corrected WIMP relic density from 
the loop-corrected WIMP annihilation cross section. 

The evolution of the WIMP number density with time t in the early Universe is 
governed by the Boltzmann equation [13] 



Here H is the Hubble parameter describing the expansion of the Universe, a is the total 
WIMP annihilation cross section, v is again the relative velocity between the annihilating 
WIMPs in their center of mass frame, n^,eq is the WIMP number density in thermal 
equilibrium, and (...) denotes thermal averaging. For non-relativistic kinematics, the latter 




(14) 
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{av) = ^ {avfe-^^y^dv, (15) 



is given by 



where x = m^/T, T being the temperature of the thermal bath. 

As long as the WIMP annihilation (or creation) rate is larger than the Hubble expan- 
sion rate, the WIMPs are (almost) in thermal equilibrium. However, once the annihilation 
rate falls below the expansion rate, WIMPs nearly decouple. Their present relic density is 
then to very good approximation given by [13] 



^ ,2 8.5 X 10-11 GeV-2 

il^h = ^ . (16) 

g*{xF)J{xF) 

Here xp = i^xf^F^ where Tp is the freeze-out temperature of the WIMPs, is the number 
of relativistic degrees of freedom, and the annihilation integral J{xf) is defined as [I3] 



J{xf) = / dx^-^ . (17) 

The freeze-out temperature of typical WIMPs is rather small, Tp ^ hence 
WIMPs are non-relativistic when they freeze out. This suggests an expansion of aov in 
powers V |13] : 

aov^£/ + Mv^ + ■ ■ ■ , (18) 

where ^ and ^ are independent of v. Note that contains only S*— wave contributions, 
while ^ contains both 5*— and P— wave contributions. While this often [9] works rather 
well for the tree-level cross section (Tq, the loop corrections we computed in the previous 
Section cannot be parameterized in this way. We saw in Fig. H] that the correction factors 
Is and Ip depend strongly on v via the quantity i/r = fi/lpl = 2fj,/{m^v). We thus have 
to re-compute (crv), which can then be used in Eq. fll7l) to derive the Dark Matter relic 
density using Eq.(IT6l). 

The 1-loop corrected WIMP annihilation cross section for partial wave labeled by L 
can be written as 

o-L = CTo.L + So-p . (19) 
Eq.dH]) and the analogous expression for annihilation from the P— wave imply that 

Sap = j^Ip{r)ao,p. (20) 

The thermal average over the tree-level cross section, expanded according to Eq.( |T8|) . can 
be computed easily [T3] : 

(aot;)(x) ~^+ — . (21) 

X 

In order to calculate the thermal averages over the correction terms, we rewrite the integral 
in Eq. ffTSl) in terms of the integration variable t = V\/x: 

{5aLv) = —^ dtt\6aLv)e-'/\ (22) 
10 



Eqs.f l20|) imply that the t— dependence of Saiv takes the form v^Iiir) = [t/ ^y''lL ^^j^ j ; 

with n = — 1 (+1) for annihilation from an S* — (P— )wave initial state. This implies that 
the thermal average over the correction terms to the WIMP annihilation cross section can 
be written clS X "'^ times a function of 2; = ly^^fx j m^. 

In Appendix A we describe how we parameterize the resulting functions for S— and 
P— wave annihilation. The one-loop corrections to the "annihilation integrals" can then 
be written as 



6Jp{xf) 




agz"^ + hsz + cs 
°° dz 



Zp Z 



exp(— ap2; + hp] 



cpz + di 



(23) 



Here zp = {2fi/m^)y/xF, and the coefficients as, bs, cs, ds and ap, bp, cp, dp are given 
in Eqs. fl3T]) and f p5|l . respectively. We have made the simplifying assumption that is 
dominated by P— wave contributions. This is true whenever the contribution oc ^ to the 
annihilation integral is comparable to, or dominates over, the one oc In the opposite 
case the correction to the e^^— term will in any case be insignificant]^ 

Not surprisingly, the corrections are quadratic in the coupling of the exchanged boson 
to the WIMP. Moreover, they depend on the boson mass only through the ratio ^i/m-^. 
Recall finally that Is has to be multiplied with —3 if the exchanged boson has axial vector 
couplings; the same factor has to be included in 6Js{xf). 

In Fig. |5] we show the relative size of the loop corrections to the annihilation integrals 
for S— and P— wave annihilation, divided by the square of the WIMP-boson coupling 
g. Eq.( IT6|) shows that this also gives the relative change of the relic density due to our 
loop corrections, as long as the corrections are small, in which case xp remains essentially 
unchanged by these corrections; note that xp depends only logarithmically on the annihi- 
lation cross section. This is the main model-independent result of our paper. The size of 
the corrections to the relic density can be read off directly, by simply inserting the values 
of masses /i and and coupling g given in a concrete model. 

We see that the corrections are less important for P— wave annihilation. This is true 
even for fi <^ m^, where the loop functions Is and Ip become equal. In this limit v6as oc 1/v 
while v6ap oc v. Performing the integrals for the thermal averaging, Eq. (fT5|) . and inserting 
the results into the definition of the annihilation integral, Eq.f ll7p . one finds 

oJs/Js 9 

The ratio becomes even smaller for nonvanishing /i, because then Ip < Is- 

Fig. [5] shows results for xp = 25 (upper curves) and xp = 20 (lower curves); this 
spans the range of decoupling temperatures in usual WIMP models. One can easily show 



^The one-loop correction to the 0{v'^) contribution to the 5— wave annihilation cross section involves 
the sum of two terms: the product of the 0{v°) one-loop correction and the 0{v'^) tree-level amplitude, 
and the product of the 0{v'^) one-loop correction and the 0{v^) tree-level amplitude. Only the first of 
these terms can be computed using the results of Sec. 2. 
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Figure 5: Corrections to the annihilation integrals for annihilation from the S— wave (solid, 
black) and P— wave (dashed, red or grey), normalized to the tree-level results, and divided 
by the square of the WIMP coupling to the exchanged boson g. The upper (lower) curves 
are for inverse scaled freeze-out temperature xp = m^/Tp = 25 (20). 

analytically that in the limit yU — 0, 6Jl/Jl oc ^Jxf for both L = (5— wave) [11] 
and L = 1 (P— wave). On the other hand, Fig. [5] shows that the relative correction to 
the annihilation integral becomes independent of xp once /i > 0.3m-;,.. This figure also 
shows that for the most plausible scenarios with electroweak strength couplings, possibly 
suppressed by mixing effects, where < 0.5, the loop corrections are significant only for 
/i < m^, as stated in the beginning of Sec. 2. 

Two comments are in order before concluding this Section. First, Fig. [5] seems to 
imply that, especially for P— wave annihilation, the corrections are never very large. This 
is misleading. A simple one-loop calculation can be trusted only \i ba < o for all relevant 
velocities v. This requires g'^Is,p / {v'n'^) < 1. For /i — )■ this will always be violated at 
sufficiently small requiring summation of higher orders. In the standard treatment jHO 
IHj this leads to 5(T oc 1/u even for annihilation from a P— wave initial state. We saw at the 
end of Sec. 2 that for finite /i the maximal size of the correction to the cross section is of 
relative size g'^m^/ {Ann). Our one-loop calculation can be trusted as long as this quantity 
is well below 1. 

Secondly, it has very recently been pointed out [U] that Eq.( lT7|) becomes inadequate 
for very small In this case the annihilation integral receives sizable contributions from 
quite large x, i.e. from low temperatures. Eq. f|T7|l assumes that the WIMPs are in kinetic 
equilibrium while they are annihilating. It has to be modified for temperatures below 
the kinetic decoupling temperature, which is typically a few (tens of) MeV [16]. This 
modification can have sizable effects for very small [I5]. However, we just saw that 
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our strictly perturbative treatment is not reliable in this case anyway. Recall that v6a 
becomes constant, rather than scaling like 1/v, for velocities below fcrit = f^/f^x- ^^^^ 
perturbative treatment will be reliable only if ii/m^ 3> g^/{4:7T) ~ 0.01. This implies that 
WIMP annihilation will quickly become irrelevant for x > l/t'crit ~ 10^5 ^^11 before kinetic 
decoupling occurs. 

4 Applications 

In this Section we apply our results to existing WIMP models. We start with two simple 
models with scalar or fermionic SM singlets; in the third Subsection we discuss the more 
widely studied case of the MSSM neutralino. 

4.1 Scalar singlet WIMP 

This is the simplest WIMP model ^7\. One only needs to introduce a single real scalar 
field X to describe Dark Matter. If one forbids terms linear in x by some (possibly discrete) 
symmetry in order to prevent x decays, the only renormalizable coupling to SM fields 
allowed by all symmetries is of the form x^l^P? where h is the scalar Higgs doublet. Upon 
weak symmetry breaking this gives rise to a trilinear scalar interaction of the form Vx'^4', 
where is the physical Higgs scalar of the SM and V = 246 GeV the vacuum expectation 
value (vev) of the Higgs. 

These interactions allow x to annihilate via exchange in the s— channel; annihilation 
into two (f) bosons is also allowed for > m^. An accurate tree-level calculation of the 
resulting relic density has been performed in [18]. Writing the coefficient of the x^l^P 
term in the Lagrangian as —k/2, they find that the correct relic density ([1]) is obtained 
for k ~ 0.28m^/(l TeV), unless ~ in which case an even smaller k is required. 

We can use our formalism to compute corrections to this result from exchange prior to 
annihilation, i.e. </) = in this case. This gives a coupling factoi0 k'^V'^ / {Am^ ~ 0.0012, 
(almost) independent of m^. Fig. [5] shows that the corrections due to exchange in the 
initial state will then at best be at the permille level. In this model we therefore do not 
find any significant radiative corrections involving the initial state only. This should also 
hold for the inert doublet model [H], since the coupling between the inert WIMP doublet 
and the SM Higgs boson will have to satisfy a similar relation. 

4.2 Fermionic singlet WIMP 

The next simplest WIMP model [20] contains a Dirac fermion SM singlet x &s well as a 
real scalar singlet y?, with couplings gxX'P + ^V'l^P) where h is again the Higgs doublet of 
the SM. The latter term induces mixing between the singlet ip and the SM Higgs boson, 
allowing ip to decay. If this mixing is small and m^p < m^, the dominant XX annihilation 

^"^Rccall from our discussion near the end of Sec. 2 that the relevant quantity for a purely scalar theory 
is the trilinear scalar coupling divided by 2-171^, i.e. in the case at hand, g = kV/(2m^). 
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channel is into two (p bosons, via t— and n— channel diagrams. In the non-relativistic limit 
this is a pure P— wave process, with tree-level cross section 



vaixx = ; ^-h + Oi^") > (25) 



where (3 = y 1 — m'^/m'^. In the limit ^ this simplifies to 



1287rm2 



The tree-level calculation therefore predicts the correct relic density ([T]) for coupling 



n m^ 



This can then be used to read off the correction due to (p exchange in the initial 
state from Fig. [5l with fi = m^p. The result is shown in Fig. |6l where we have used the 
exact (tree-level) result ( 125|) to derive the required coupling strength. The black curves 
show that this coupling strength depends only very weakly on the mass of the scalar for 
fi = m^ < 0.5m^. Note that the cross section slightly increases with increasing m^ as 
long as < 0.85m^. The reason is that increasing allows the t— and m— channel 
propagators to be less off-shell, as shown by the denominator in Eg. (1251) . As a result, the 
required coupling strength slightly decreases with increasing m^. 

However, the red curves show that this effect is much smaller than the strong de- 
pendence of SJp/Jp on fJ'/m^ illustrated in Fig. [51 More importantly, we see that the 
corrections due to (f exchange in the initial state can easily exceed the uncertainty of the 
observational determination of the DM relic density; for fJ^/m^ = 0.5 (0.2, 0.1) they even 
reach the 10% level for m^ > 450 (210, 130) GeV. Since the correction is positive, one 
would have to reduce the coupling in order to obtain the correct relic density after inclusion 
of one-loop corrections. This would correspondingly reduce all interactions between the 
WIMP X and the SM particles, all of which are mediated by ip exchange. 



4.3 The lightest neutrahno in the MSSM 

We finally want to apply our formalism to the lightest neutralino in the minimal supersym- 
metric extension of the Standard Model (MSSM), which is the probably best motivated 
WIMP, and certainly the most widely studied [H [21] one. For simplicity we will assume 
that sfermions are heavy. Given experimental lower bounds on the masses of sfermions and 
Higgs bosons, relatively light sfermions by themselves typically only lead to an acceptable 
neutralino relic density in the presence of significant co-annihilation [22]. This involves 
several particles, with mass splittings of order of the absolute value of the 3-momentum in 
the initial state. These more complicated scenarios cannot be treated with the formalism 
presented in this paper0 

^^Very recently the summation of "Sommerfeld corrections" for the case of nearly degenerate states was 
discussed in [23]. 
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Figure 6: Strength of the XXV^ couphng g (black), as well as the resulting one-loop correc- 
tion to the annihilation integral (red), in a model with Dirac fermion singlet DM. g has 
been computed from the requirement that XX ~^ V^V^ annihilation produces the correct x 
relic density at tree level. The sohd (dashed, dotted) curves are for n/m^ = 0.5 (0.2, 0.1), 
where fi = m^. 

Generally speaking, in the MSSM the neutralinos are mixtures of the U{1)y gaugino 
B, the neutral SU{2) gaugino W^, and of the two neutral higgsinos h^, /ig- 

X° = NaB + N,2m + N,ih', + N^Ji^ (z = 1, ■ ■ ■ , 4). (26) 

The coefficients Ni^ satisfy the sum rule Ylt=i \^ik\^ = 1 Most phenomenological analy- 
ses of the MSSM assume that the soft SUSY breaking gaugino masses unify at or near the 
scale of Grand Unification [21]. This implies that the f/(l)y gaugino mass is about half the 
SU{2) gaugino mass near the TeV scale. As a result, the Wino component of our candidate 
WIMP, the lightest neutralino (x = Xi); is subdominant, i.e. jA?"!!!^ ^ 1-^121^. If sfermions 
are heavy, x annihilation involves couplings of the lightest neutralino to gauge or Higgs 
bosons, which vanish in the pure Bino limit (|A^ii| — J- 1). In models with gaugino mass 
unification and heavy sfermions, the annihilation cross section can thus only be sufficiently 
large if x has significant higgsino components. 
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On the other hand, for a nearly pure higgsino, where |A'^ii| + |iVi2| ^ 1, the cross 
section for annihilation into W^W~ and Z^Z^ pairs is so large that a mass ^ 1 
TeV is required to obtain the correct relic density ([1]) |23]. Such a large mass for the 
lightest superparticle is at odds with the primary motivation for postulating the existence 
of superparticles, which is the stabilization of the weak scale against quadratically divergent 
quantum corrections. 

Assuming gaugino mass unification, the most natural neutralino satisfying the con- 
straint ([T]) is therefore a bino-higgsino mixture, dubbed a "well-tempered neutralino" in 
ref.[2S]. Note that the neutralino couplings to Higgs bosons involves products of a combi- 
nation of gaugino components {Ni2 — tan^vy-^ii, where 9w is the weak mixing angle) with 
one of the higgsino components (A^is, No). These couplings are maximal in the region 
of strong gaugino-higgsino mixing, and should thus be sizable for the "well-tempered" 
neutralino. Moreover, as well known, at least one of the neutral MSSM Higgs bosons is 
rather light, with mass below 130 GeV. This leads one to expect potentially sizable 1-loop 
corrections due to the exchange of a Higgs boson prior to WIMP annihilation. 

We checked this with the help of the the code micrOMEGAs 2.2 [SB]. Among other 
things, this program computes the complete tree-level neutralino annihilation cross sec- 
tions for all two-body final states. It does not resort to the non-relativistic expansion f lTSj) . 
but one can easily determine the coefficients and ^ by calculating the annihilation cross 
section at two different values of i.e. for two (slightly) different cms energies i/i. These 
coefficients are then used in Eqs. ( l23l) to compute the corrections to the annihilation inte- 
gral; we emphasize that we continue to use the full cross sections, not their non-relativistic 
expansions, for the calculation of the tree-level contribution to the annihilation integral. 
We also take xp from the program. The one-loop corrected x I'^lic density can then be 
expressed as 

n^h-" = fi^,o/,2^M^£L^ (27) 

Jo{xf) + 6J{xf) 

where 6 J = 6Js+SJp, and the tree-level value Jo{xf) can be calculated from the program's 
tree-level prediction fixfi^^ using Eq. (fTB]) . 

The result is shown in Fig. [3 It has been generated using micrOMEGAs [26], using 
Softsusy [27] to calculate the superparticle and Higgs boson spectrum; we specify the input 
directly at the weak scale. All points satisfy the relic density constraint ([T]) within two 
standard deviations. For simplicity we take (unnaturally) large values for the masses of all 
sfermions and of the CP-odd Higgs boson A, but the results would not change significantly 
as long as rrigf > 2m^ and > 3m^. We include loop corrections due to exchange of the 
Z— boson as well as both CP-even neutral Higgs boson, but the contribution from the 
heavier Higgs boson is totally negligible due to its large mass. 

We see that the corrections are most important for near 350 GeV. For smaller 
WIMP mass the corrections are reduced because the higgsino component of x becomes 
smaller, and because the ratio of light Higgs and WIMP masses becomes smaller, which 
reduces the loop functions. The latter effect would tend to increase the correction for 
heavier WIMPs. However, for > 350 GeV the gaugino components of x decrease 
quickly; this reduces its coupling to the light Higgs boson, which is most important here. 
Moreover, for > 900 GeV, co-annihilation with xt X2 become important [TH |2S1 
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tanp=2.5, msf=5TeV, mA=5TeV 




Figure 7: The relative size of one-loop correction to the relic density of a "well-tempered" 
neutralino. The ratio between Mi and /i is fixed such that the relic density satisfies the 
constraint ([1]) within two standard deviations. We take M2 = 2Mi, consistent with gaugino 
mass unification, tan/3 = 2.5, and assume sfermions and most Higgs bosons to be very 
heavy. Scenarios with < 200 GeV are excluded by Higgs boson searches. 

|24] . The effect of light boson exchange corrections to co-annihilation is beyond the scope 
of our paper, and has not been included in Fig. [71 As a result, the correction becomes 
comparable to the anticipated post-PLANCK precision of the observational determination 
of Q^h"^ only for a rather narrow range of ~ 350 GeV. This is consistent with the results 
of Fig. |5l given the fact that the coupling of our "well-tempered" neutralino to the lightest 
Higgs boson does not exceed 0.2. 

The small size of the corrections due to boson exchange in the initial state indicate 
that these may well not be the leading radiative corrections in the MSSM. In fact, full 
electroweak one-loop calculations [TH [22] found much larger corrections in some cases. 
These are presumably due to UV-sensitive effects, which cannot be treated using our 
formalism. Moreover, QCD corrections can significantly affect the annihilation cross section 
into quarks [TU [30] . 
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5 Summary and Discussion 



We have calculated one-loop corrections to the WIMP annihilation cross section, and the 
corresponding corrections to the relic density, due to the exchange of a relatively light 
boson between the WIMPs. 

The formalism for calculating the corrections to the cross section is described in Sec. 2. 
We used a non-relativist ic formalism [TOl |5], i.e. we only included small loop momenta. 
The correction can then also be understood as re-scattering of the WIMPs prior to their 
annihilation. The motivation for this is that these configurations can give rise to sizable 
corrections if the exchanged boson is lighter than the WIMP and the relevant coupling is 
not too small. Note that these corrections are universal, i.e. they are the same for all final 
states. However, we saw that they do depend on the partial wave of the initial state, being 
smaller for annihilation from the P— wave if the exchanged boson is massive. The main 
result of this Section is that these corrections can be described by loop functions which 
only depend on the ratio of the mass of the exchanged boson and the cms three-momentum 
of the annihilating WIMPs; simple yet accurate parameterizations for these loop functions 
are given in Eqs.f llSp . We also checked explicitly that the correction is independent of the 
form of the coupling of the exchanged boson, except for the case of axial vector exchange 
in the S*— wave, where an additional factor of —3 is required. Moreover, the corrections can 
be used for scalar as well as fermionic WIMPs. 

In Sec. 3 we computed the resulting correction to the relic density. We saw that the 
thermal averaging and the calculation of the "annihilation integral" further reduce the 
relative importance of the corrections in case of annihilation from a P— wave initial state. 
The main model-independent result of our paper is shown in Fig. |5l which shows the 
relative size of the corrections to the annihilation integral - or, almost equivalently, to 
the relic density - as a function of the ratio of the masses of the exchanged boson and 
the WIMP. We saw that for 0{1) coupling between boson and WIMP, the corrections are 
significant unless the boson is heavier than the WIMP. 

Finally, in Sec. 4 we applied this formalism to several WIMP models. We showed 
that these corrections are very small for a scalar singlet WIMP, but can be very large for a 
Dirac fermion singlet WIMP annihilating into a light scalar singlet. Finally, we analyzed the 
"well-tempered" neutralino, which is a mixture of U{1)y gaugino and higgsinos, and found 
corrections comparable to the anticipated post-PLANCK precision of the observationally 
determined WIMP relic density only for a narrow range of neutralino masses near 350 
GeV. The main reason is that the relevant coupling is always below 0.2 in this case. In a 
supersymmetric scenario, bigger couplings of WIMPs to Higgs bosons are e.g. possible for 
"singlino" Dark Matter in scenarios where the MSSM is extended by an additional Higgs 
singlet superfield [3Tj . 

In the MSSM diagonal couplings of the lightest neutralino to gauge or Higgs bosons are 
always suppressed by mixing angles. Off-diagonal couplings may be large, however. These 
can give rise to sizable corrections if these couplings involve particles that are only slightly 
heavier than the lightest neutralino. Such states, with sizable couplings to the lightest 
neutralino, often exist in regions of parameter space where co-annihilation is important. 
In order to treat this, one has to extend the formalism presented here to scenarios where the 
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particles in the loop have slightly different masses than the annihilating external particles. 
Moreover, if co-annihilation with a sfermion is important, there are vertex corrections 
involving the exchange of an SM fermion, rather than a boson, between the co-annihilating 
particles. This offers another avenue for future work. 

Acknowledgment 

We thank A. Pukhov for his help with using micrOMEGAs. KIN is supported in part 
by Grants-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, 
Science, and Technology of Japan, and by the Grant-in-Aid for Nagoya University Global 
COE Program, "Quest for Fundamental Principles in the Universe: from Particles to the 
Solar System and the Cosmos", from the Ministry of Education, Culture, Sports, Science 
and Technology of Japan. MD and JMK are partially supported by the Marie Curie Train- 
ing Research Networks "UniverseNet" under contract no. MRTN-CT-2006-035863, and 
"UniLHC" under contract no. PITN-GA-2009-237920. 



Appendix A: Parameterizations of the thermal aver- 
ages 



The thermal average of the correction to the 5"— wave annihilation cross section is given by 



has been defined via the nonrelativistic expansion of ov in Eq.( IT8|) . The f— dependence 
of Is can be read off Eq. (fT3!) : 



Here we have introduced the quantity u = 2fi/m^. We showed in Sec. 3 that the integral 
in Eq.f l28p is a function of the variable z = u^fx. We find the following fitting function for 
the "thermally averaged" {Is/v), defined as the integral in the last line of Eq.f l28p : 




(28) 




(29) 




(30) 



with 



as = 0.000593; bg = 0.03417; cs = 0.1015; ds = 0.1182. 



(31) 



For the P— wave. 
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^ 1^ g-x.V4^ (32) 



4vr5/2 
with 

^ ^5^' ^ - - . (33) 

1 I 3u , ^2 , > 2.1 

As above, we find a fitting function for tlie "tliermally averaged" (f/p): 

{vipht = ^ (e--^-^- + ^—-) , (34) 

\ cpz + dp/ 

with 

ap = 0.318; 6p = 0.1226; cp = 0.3309; dp = 0.6306. (35) 



Appendix B: Comparison to a full one— loop calculation 

In this Appendix we compare our approximate treatment of corrections due to ip exchange 
with a full one-loop calculation. We do this in the framework of a purely scalar theory, 
where the exact vertex correction is UV finite. As we remarked in Sec. 2, our formalism 
will not capture corrections associated with the renormalization of the coupling(s) relevant 
for WIMP annihilation, so chosing an example with UV-finite vertex correction greatly 
simplifies the comparison to the full one-loop calculation. 



'^l • • • 

; A /"'"'< V 

,V / 

x/ x/ x/ 

a) b) c) 



Figure 8: Feynman diagrams describing full initial state radiative corrections: vertex correction 
(a), wave function renormalization (b) and real emission (c); wave function renormalization of, 
and emission off, the upper leg also has to be included. The WIMP x and the exchanged boson 
are denoted by black and red (grey) dashed lines, respectively, while the blob denotes the XX 
annihilation vertex, which is independent of the x momenta. 



The Feynman diagrams describing exact one-loop corrections associated with the initial state 
are shown in Fig. 8. Here the blob describes the (tree-level) x annihilation process; this could e.g. 
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be a quartic vertex involving two lighter scalars, or a trilinear vertex coupling to the s— channel 
propagator of another scalar particle. For the purpose of our calculation we only need to know 
that the rest of the diagram described by the blob is independent of the loop momentum. We 
describe the ^XX vertex by the (dimensionful) coupling k. 

Let us begin by computing the vertex correction; recall that this is the only diagram that 
contributes in the approximate treatment of Sec. 2. It gives: 

^vertex ■ 2 /" 111 (rip^ 

"AT ~ i (2^)4 (P + qY -ml{P- qY - m2 (p - g)2 - ^2 ' ^ > 

Here is the tree-level matrix element described by the blob in Fig. 8. Recall that P = 
{Pi +P2)/2, p = {pi — P2)/2, where pi^2 are the 4-momenta of the incoming WIMPs, // is the 
mass of and k is the XX^ coupling. 

The loop integral in Eq.(|36p can be computed straightforwardly using Feynman parameters, 
giving 

^vertex ^ ri f 2 2\ /nnx 

Aq =~Y^'-^o{s,m^,m^,m^,fi). (37) 

Here Cq is the scalar Passarino-Veltman three-point function in the convention of ref. |32j . 

The loop integral in Eq. (f36l) can also be evaluated directly, following the steps of Sec. 2 



but without making any approximations in the propagators. We first perform the energy (go) 
integrals by contour integration, by summing over the residues of all poles in the lower half 
plane. In general, there are three such poles: 

pole 1 x-) 

% = w - Po ; 



pole 2 I 7-, 



C = ^J{p-q)' + ^^', (38) 

where uj = \J(f' + 7ti| as in Eq.([5]). Only the first pole has a residue that diverges in the limit 
i?, g — 7- 0. The third pole comes from the energy dependence of the (p propagator, which has been 
ignored in the approximate treatment of Sec. 2. The angular integrals can also be performed 
straightforwardly. After some algebra, we arrive at: 



^vertex ^ 1 



^0 IQtt"^ APq\p\ J Q 



\mq\ 



i2 



u:{u - Po) {\p\ - \c^f + /.2 - (a; - P,f 
1 , {\p\ + \q\f + ^^'-{co + P,^' 



u;{u: + Po) (|p-l_|g-l)2 + ^2_(^ + pg)2 
1 f^^ /i 2 + 2\p\\q] - 2Pou^ + 2\p\\q\ + 2Pou;A 

2-2|pl|g-l -2PoOJ^ ''/i2-2|pl|g-1 + 2Po^^/ 




In the last line, we have introduced uj^p = \/ q^ + ^j?. 

It is easy to see that the first term reduces to our expression ([7]) if we use the non-relativistic 
expansion for tj, which includes dropping the terms [oj — Pq)"^ in the argument of the logarithm. 
However, for large |g| these latter terms are important. They imply that the logarithm approaches 
the constant value 2|p|/Po in the limit |g| — )• oo, rather than vanishing as in Eq.([7]). As a result, 
the first line of the right-hand side (rhs) of Eq. ()39p by itself is logarithmically UV-divergent. The 
second line contributes the same UV divergence again; only after adding the contribution in the 
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third line we obtain a UV-finite result. This third line comes from the third pole in Eq. (l38p . which 
does not exist if one drops the energy dependence of the (p propagator, as in Sec. 2. This proves 
our statement in Sec. 2 that omission of the energy dependence of the ip propagator necessitates 
the use of a non-relativistic expansion in the argument of the loop integral. 

While the third line in Eq. ()39p is necessary to obtain a UV-finite result, it introduces a new 



problem: for /x — )■ it becomes IR-divergent! In the limit /i^ 
Eq. ([39]) simplifies to 

l6Po|plk1 1 



^ mr,, \p\ <^ Pq the last line of 



k1. 



As noted earlier, the UV 



The d\(f\ integration will then lead to a negative term oc In 
divergence for l^max oo precisely cancels those from the first two terms in Eq. (|39p . The 
resulting term oc In ^ in the exact vertex correction is IR-divergent for /i — t- 0. This term 
becomes significant for small fi, especially if the velocity v is not too small. This explains our 
statement in Sec. 2 that our approximation does not describe the exact vertex correction very 
well for small ^. 

However, the IR divergence does not exist in the full one-loop calculation. We have to add 
wave function renormalization (Fig. 8b) as well as real emission diagrams (Fig. 8c) in order to 
obtain an IR-finite result for /x — )• 0. This implies that adding these additional contributions 
should also remove all terms oc In ^ in the complete one-loop corrected cross section. Using 
on-shell renormalization for the wave function renormalization constant is finite!^ 



IGvr- 



(40) 



where Bq is the derivative of the scalar Passarino-Veltman two-point function with respect to 
its first argument. Adding this negative contribution doubles the IR-divergence for /i — )• 0. 

Finally, we have to treat the real ip emission diagram of Fig. 8c, plus the contribution where 
ip is emitted off the other x line. Writing the 4-momentum of the emitted ip as k = {ko, k), we 
have 



I ^real 



\An 



2(//2_2A;oPoj 



(/i2-2A:oPo)'-4(fc-p)^ 



(41) 



Performing the angular integrations of the (p phase space, this gives 



f^real 



K 

2^ 



dkn 



+ 



{2Poko - 

1 



2^2 



4/,2p2 

2A:oPo 



In- 



fi^ + 2\k\\p\ 



4\p\ {2Poko - /i2) 2k.P 



2\k\\p\ 



(42) 



where 



^2 and /cq,] 



Pq + /i^/(4Po)- In the limit /x — this also produces a 



logarithmic IR divergence, this time with positive sign, which (not surprisingly) cancels the sum 
of the IR divergent terms from the vertex correction and wave function renormalization. 

Fig. [9] shows that the sum of the vertex correction, wave function renormalization and real 
emission contributions very closely matches our approximate result of Sec. 2 for jj, < 0.5m^. In 
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'Diagram 8b also gives a logarithmically divergent contribution to ra^. This is simply removed by the 
mass counterterm in on-shell renormalization. 
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Figure 9: Normalized higher order contributions to the x annihilation cross section for rela- 
tive initial state velocity v = 0.1. The solid (black) curve shows our approximation of Sec. 2, 
given by Is /v. The dashed (blue), dotted (red) and dot-dashed (green) curves show exact 
contributions from vertex corrections, real emission diagrams, and wave function renor- 
malization, respectively; the latter has to be multiplied with —1. The dash-doubledotted 
(magenta) curve shows the sum of these three contribution, i.e. the exact one-loop correc- 
tion; it nearly coincides with the black curve for yu/m^ < 0.2. 

fact, the difference is always of order /{^/(Svr^m^), without any potentially large factors like 1/v 
or ln(m;^/ ^). We checked that this remains true at least for all v < 0.5 relevant for the calculation 
of the relic densitvl^ It is not surprising that our approximate treatment does not treat such 
"generic" higher order contributions correctly. However, our approximation does closely resemble 
the exact result whenever the latter is large. This is all we aspired to, and in most cases all we 
need when calculating DM relic densities, even if PLANCK data reduce the uncertainty of the 
observed value to the percent level. 

^■^Our approximation might fail badly when the annihilating WIMPs become ultra-relativistic, but this 
is of no concern in the present context. 
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